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ABSTRACT 

We discuss the 21cm power spectrum (PS) following the completion of reionization. 
In contrast to the reionization era, this PS is proportional to the PS of mass density 
fluctuations, with only a small modulation due to fluctuations in the ionization field 
on scales larger than the mean-free-path of ionizing photons. We derive the form 
of this modulation, and demonstrate that its effect on the 21cm PS will be smaller 
than 1% for physically plausible models of damped Lyev systems. In contrast to the 
21cm PS observed prior to reionization, in which HII regions dominate the ionization 
structure, the simplicity of the 21cm PS after reionization will enhance its utility 
as a cosmological probe by removing the need to separate the PS into physical and 
astrophysical components. As a demonstration, we consider the Alcock-Paczynski test 
and show that the next generation of low- frequency arrays could measure the angular 
distortion of the PS at the percent level for z ~ 3 — 5. 

Key words: cosmology: diffuse radiation, large scale structure, theory - galaxies: 
high redshift, intergalactic medium 



1 INTRODUCTION 

Recently, there has been much interest in the feasibility of 
mapping the three-dimensional distribution of cosmic hy- 
drogen through its resonant spin-flip transition at a rest- 
frame wavelength of 21cm (Furlanetto, Oh & Briggs 2007; 
Barkana & Loeb 2007). Several experiments are currently 
being constructed (including MWA Q, LOFAR 0, PAPER 
0, 21CMA0, GMRT0) and more ambitious designs are be- 
ing planned (SKA Ej- 

One driver for mapping the 21cm emission is the possi- 
bility of measuring cosmological parameters from the shape 
of the underlying power spectrum (PS; see Loeb & Wyithe 
2008). During the epoch of reionization, the PS of 21cm 
brightness fluctuations is shaped mainly by the topology of 
ionized regions, rather than by the PS of matter density fluc- 
tuations which is the quantity of cosmological interest (Mc- 
Quinn et al. 2006; Santos et al. 2007; Iliev et al. 2007). As 
a result, the line-of-sight anisotropy of the 21cm PS due 
to peculiar velocities must be used to separate measure- 
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ments of the density PS from the unknown details of the 
astrophysics (Barkana & Loeb 2005; McQuinn et al. 2006). 
The situation is expected to be simpler both prior to the 
formation of the first galaxies (at redshifts z > 20, Loeb 
& Zaldarriaga 2004; Lewis & Challinor 2007; Pritchard & 
Loeb 2008), and following reionization of the intergalactic 
medium (IGM; 1 < z < 6) - when only dense pockets 
of self-shielded hydrogen, such as damped Lya absorbers 
(DLA) and Lyman-limit systems (LLS) survive (Wyithe & 
Loeb 2008; Chang et al. 2007; Pritchard & Loeb 2008). 

In this paper we focus our discussion on the post- 
reionization epoch (Wyithe & Loeb 2007; Chang et al. 2007). 
The DLAs which contain most of the neutral hydrogen 
mass in the Universe at z < 6 are expected to be hosted 
by galactic mass dark matter halos (Wolfe, Gawiser & 
Prochaska 2005). A survey of 21cm intensity fluctuations 
after reionization would measure the modulation of the 
cumulative 21cm emission from a large number of galax- 
ies (Wyithe & Loeb 2008; Wyithe, Loeb & Geil 2008; Chang 
et al. 2007). Regarding the measurement of the 21cm PS, 
this lack of identification of individual galaxies is an advan- 
tage, since by not imposing a minimum threshold for de- 
tection, such a survey collects all the available signal. This 
point is discussed in Pen et al. (2008) , where the technique is 
also demonstrated via measurement of the cross-correlation 
of galaxies with unresolved 21cm emission in the local Uni- 
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Studying the 21cm PS after (rather than during) reion- 
ization offers two advantages. First, it is less contaminated 
by the Galactic synchrotron foreground, whose brightness 
temperature scales with redshift as (1 + z) 2 ' 6 (Furlanetto, 
Oh & Briggs 2006). Second, because the UV radiation field is 
nearly uniform after reionization, it should not imprint any 
large-scale features on the 21cm PS that would mimic the 
cosmological signatures. In addition, on large spatial scales 
the 21cm sources are expected to have a linear bias analo- 
gous to that inferred from galaxy redshift surveys. 

Most previous studies of post-reionization 21cm fluctu- 
ations have assumed that the 21cm emission traces pertur- 
bations in the matter density, and have considered peculiar 
motions only after a spherical average (Wyithe & Loeb 2008; 
Pritchard & Loeb 2008, Loeb & Wyithe 2008). The excep- 
tion is a recent paper discussing measurements of the dark 
energy equation of state (Chang et al. 2007). Since the neu- 
tral gas resides within collapsed dark matter halos, galaxy 
bias plays an important role in setting the 21cm fluctua- 
tion amplitude. In this paper we derive the 21cm PS af- 
ter reionization in the context of the formalism that has 
been developed to calculate it during reionization (Barkana 
& Loeb 2005; McQuinn et al. 2006; Mao et al. 2008). By 
framing the derivation this way, the relative merits of cos- 
mological constraints from 21cm surveys at redshifts before 
and after reionization can be more easily understood. In ad- 
dition, this formalism provides a framework to describe the 
possible effect of fluctuations in the ionizing background. 
We then compute the Alcock-Paczynski effect (Alcock & 
Paczynski 1979) as an example for the cosmological utility 
of the post-reionization 21cm PS. In our numerical exam- 
ples, we adopt the standard set of cosmological parameters 
(Komatsu et al. 2008), with values of Q b = 0.24, fi m = 0.04 
and £Iq = 0.76 for the matter, baryon, and dark energy 
fractional density respectively, and h — 0.73, for the dimen- 
sionless Hubble constant. 



2 21CM POWER SPECTRUM AFTER 
REIONIZATION 

The 21cm PS after reionization is expected to be domi- 
nated by the neutral content of galaxies. In a scenario where 
fluctuations in the ionizing background can be ignored, the 
form of the PS could therefore be inferred directly from the 
galaxy PS (Loeb & Wyithe 2008). However, we also include 
here the possible influence of a fluctuating ionizing back- 
ground on the 21cm PS. Throughout our discussion, we as- 
sume the baryonic overdensity 5 to equal the overdensity in 
dark-matter on sufficiently large scales. 



2.1 The 21cm brightness temperature 

The 21cm brightness temperature fluctuation in a region of 
IGM is 

AT = 23.8 ( [1 - Xi(l + S x )] (1 + 5) (1 - 5 V ) mK, 

(1) 

where 5 V = dv r /dr(Ha) 1 and dv r /dr is the gradient of the 
peculiar velocity along the line-of-sight. The quantities Xi 
and S x are the mean ionization fraction and the fluctuation 



in ionization fraction, respectively. We have assumed that 
the spin temperature of hydrogen is much higher than the 
CMB temperature [as expected for collisionally coupled gas 
in collapsed objects, and observed in some DLAs (Curran et 
al. 2007)]. Thus we may neglect fluctuations in the kinetic 
and spin temperatures of the neutral hydrogen gas through- 
out the post reionization epoch. To evaluate the brightness 
temperature fluctuation we need to compute 5 X . With this 
goal in mind, we first write the fluctuation in neutral hydro- 
gen fraction as 



5m = xm/xm — 1, 



(2) 



where xm and xm are the cosmic mean and local values for 
the neutral fraction (mass averaged). 

2.2 Effect of the ionizing background 

Next we compute the effect of the ionizing background on 
the neutral fraction. Two regimes must be considered. 



2.2.1 Systems which are optically-thin to ionizing 
radiation 

Low density regions of the Universe contain optically thin, 
Lye* absorbers. The ionization fraction in this regime is con- 
trolled by the balance between the ionization rate owing to 
the UV background and the recombination rate at the lo- 
cal gas density. Given an ionization rate T, and a hydrogen 
density iih with a neutral fraction xm, the equilibrium in 
the optically thin regime is given by the condition 



xm 



(3) 



where qb is the case-B recombination coefficient (Oster- 
brock & Ferland 2006). Assuming that the low density re- 
gions of the hydrogen density field are unbiased with respect 
to the underlying dark matter density, we have 



xm oc (1 + o — ), 



(4) 



where Ar is the perturbation in the intergalactic ionizing 
background. Note that the neutral fraction in optically thin 
regions is increased in overdense regions by recombinations, 
and decreased by the presence of an excess ionizing back- 
ground. 



2.2.2 Optically-thick absorbers 

Since optically-thick systems are self- shielded, the effect of 
the UV background on their ionization state depends on 
the gas distribution within them. Below we discuss the ef- 
fect of an ionizing background on the HI content of self- 
shielded systems. Since most of the hydrogen in such sys- 
tems is known to be contained within DLAs, we focus our 
discussion on DLAs. 

The nature of DLAs is not understood, although they 
are thought to be formed by dense, self-shielded gas in galax- 
ies (Wolfe, Gawiser & Prochaska 2005). The detailed mod- 
eling of DLAs is very uncertain, requiring numerical sim- 
ulations that go beyond the scope of this paper. However 
we can utilise results of simulations in the literature to esti- 
mate the range of possible strength for the influence of the 
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UV background. The column density distribution of DLAs 
has a power-law form with a feature at ~ 10 20 cm -2 (e.g. 
Storrie-Lombardi & Wolfe 2000). This feature has been in- 
terpreted by a number of authors as being due to the effects 
of self shielding of DLAs in a meta-galactic ionizing back- 
ground (e.g. Corbelli, Salpeter & Bandiera 2001; Zheng & 
Mirada-Escude 2002). 

Zheng & Mirada-Escude (2002) considered the col- 
umn density distribution of damped Lya systems modeled 
as spherical isothermal spheres. They solved for the self 
shielded HI density profile numerically and found that the 
neutral fraction reaches ~ 10~ 3 at the radius where the sys- 
tem becomes optically thick to meta-galactic ionizing radi- 
ation. Their results demonstrate that the transition from 
highly ionized to highly neutral gas occurs over a narrow 
region which is very small compared with the size of the 
absorbing system. Moreover, they find that the radius r ss 
at which the gas becomes self-shielding is proportional to 
p-V3 p or a spherical isothermal sphere, the hydrogen mass 
enclosed within radius r scales as 



M H (< r) oc r, 
and so the HI mass scales as 

M m = M H (< r SB ) oc r ss oc F 



-1/3 



(•>) 



(6) 



We therefore find that a given perturbation in the ionizing 
background produces a corresponding perturbation in the 
mass of neutral hydrogen within the self-shielded system of 
the form 



AM H 
Mm 



1 AT 

3 r ' 



(7) 



Alternatively, the distribution of gas in a DLA may 
be better represented by a self-gravitating disk. The effect 
of a meta-galactic background on the ionization structure 
and starformation has been investigated by a number of 
authors (Corbelli, Salpeter & Bandiera 2001; Schaye 2004; 
Susa 2008). The hydrogen column depth at which the gas 
becomes self-shielding was shown to be proportional to T 1 / 3 
as in the spherical isothermal case (Susa 2008). However 
the exponential profile of the disk models implies that the 
ionized portion of the gas contains much less mass than in 
the case of a more extended power-law isothermal profile. 
Corbelli, Salpeter & Bandiera (2001) compute the column 
depth of HI as a fraction of the column depth of hydrogen, 
at different intensities of ionizing background. For a disk 
with a column density of 10 21 cm -2 , the fluctuation in HI 
mass owing to a fluctuation in the ionizing background is 
approximately 



hav<Q 



xhi oc 



(i + b<5)(i-CAr/r) 

1 + 5 

l + bS-5-C 



AT 



(9) 
(10) 



where the last equality has reduced the expression to lowest 
order in 5. 



2.2.3 Combined neutral content of the IGM 

The above estimates of mass fluctuation suggest that the 
neutral hydrogen content of the IGM is modified owing to 
fluctuations in the ionizing background (Sj = AT/T) by a 
factor of the form 



(11) 



AMm _ 3 Ar 

— — — ~ -2 x 10 3 — . 
Mm T 



(8) 



where C is a constant describing the magnitude of the effect. 
For optically thin regions C = 1, while for optically thick 
absorbers C can take a range of values. In the case of an 
isothermal sphere, a fluctuation in the ionizing background 
leads to a fluctuation in the HI mass of a self shielded system 
that is of comparable magnitude, indicating that the value of 
C in equation (|11[) would be of order unity in this case. How- 
ever, in halos with virial temperatures larger than ~ 10 4 K, 
hydrogen with an isothermal profile would be collisionally 
ionized, and hence have cooled into a much more concen- 
trated profile. Since the host masses of DLAs are known 
to be greater than 10 10 Mq from clustering studies (Cook et 
al. 2006) , an isothermal profile does not provide a physically 
plausible model. On the other hand, modeling DLAs using 
an exponential disk, as is appropriate for gas rich galaxies 
implies a much smaller value value of C ~ 10 -3 -10~ 2 . Thus, 
the effect of fluctuations in the ionizing background on the 
mass averaged HI density is expected to be small (< 1%). 

We can now add the effects of fluctuations in ionizing 
background on the neutral content of optically thin and op- 
tically thick absorbers within a region of large scale over- 
density 5. We have 

xhi = Fthin5in(l + S - 5j) + .Ftjiidk*Hi(l + b5 - 5 - G8j) 

= Xm [1 + (F thi n + -Fthick(& - 1))<5 - (Fthin + CF thick )5j] 

~ xm [1 + {b - 1)5 - CSj] , (12) 

where Fthin and fthick are the fractions of HI in the opti- 
cally thin and optically thick regimes respectively. In the 
last equality, we assume the case where the optically thin 
component contains a negligible amount of the cosmic HI. 
Note that we ignore the optically thin component of the hy- 
drogen in the remainder of this paper. However we point 
out that the inclusion of this additional component changes 
only the values of the co-efficients of 5 and Sj in the above 
equations, but not the form of the expression. 



This effect is two orders of magnitude smaller than in the 
case of a spherical isothermal sphere. 

To evaluate the fluctuation in neutral fraction associ- 
ated with DLAs we assume that the neutral gas is hosted 
within halos of mass M with associated galaxy bias b. We 



7 The factor of (1 + <5) in the denominator is present because 
this is the fluctuation in neutral hydrogen fraction rather than in 
neutral hydrogen mass. 
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2.2.4 Scale dependence of fluctuations in the ionizing 
background 

Before deriving the fluctuation in neutral fraction we next 
need to calculate the dependence of the parameter Sj on 
scale. Ionizing radiation in the IGM has a finite mean-free- 
path which increases with time following the end of reioniza- 
tion. On scales larger than the mean-free-path (A m f p ), with 
associated wavenumber k < fc m f P = 27r/A m f p , all ionizing 
photons within a fluctuation were produced by galaxies that 
were also within that fluctuation. In this case fluctuations in 
the ionizing background simply trace fluctuations in the den- 
sity of galaxies so that the value of Sj oc b5 is independent of 
scale. However, on scales smaller than the mean-free-path, 
ionizing photons were not produced by galaxies local to the 
fluctuation. As a result Sj/S is scale dependent. 

To derive this dependence we convolve the real space 
density field with a filter function to account for the effects of 
finite mean-free-path and the inverse square dependence of 
ionizing flux. The fluctuation in flux at a position x becomes 



Sj(x) oc / dx'G(x,x)6(x). 



(13) 



For simplicity we assume that all ionizing photons travel one 
mean-free-path, hence 



G(x, x 



.,, _ 6(|a; - x'\ - A m fp) 



(14) 



where © is the Heaviside step function. This may be rewrit- 
ten 



5j(k) oc g(k)5k, 



(15) 



where 5j(k) and g are the Fourier transforms of the ionizing 
flux field and filter function respectively. We have 



where 



9{k) 



Si(ar) 



Si(fcA mfp ) 



A: 



/" 

Jo 



du 



sin(u) 



(16) 



(17) 



and k = \k\. For small scales fcA m f p 3> 1, we obtain Sj(k) oc 
5k/k. For large scales we find 5j(k) oc 5k as expected from 
the argument described above. 

Thus, we get the following dependence for the function 
s, which can now be written in terms of 8, 



s(S) = [1 - K(k)bS] . 



(18) 



where we have explicitly written in the galaxy bias (noting 
that it is the galaxy density rather than the matter density 
which sources the ionization field) , and defined the function 



K(k) = Ko 1 + 



mfp 



(19) 



where K oc C is a new constant, which from the discus- 
sion following equation (|11[) is expected to be smaller than 
~ 0.01. This fitting function interpolates smoothly between 
the limiting behaviour on small and large scales. We empha- 
sise that this formulation is only valid on scales for which 
fluctuations in the density field are in the linear regime. 

We can now express the fluctuation in the neutral hy- 
drogen content, noting that this depends both on the den- 



sity of galaxies and on the fluctuations in the ionizing back- 
ground. 

fat m [6(1 -*■(*)) -1] S, 

Similarly, we write the fluctuation in ionized fraction in 
terms of the local ionized hydrogen fraction (xi) as 



5 X = Xi/xi - 1. 



(20) 



Using the fact that Xi + xm ~ 1, we then find the ionization 
fluctuation in terms of the galaxy bias 

5 X « [x m /(xm - 1)] [6(1 - K{k)) - 1] S. (21) 



2.3 The 21cm PS 

We can now turn to calculating the 21cm PS. To lowest order 
in Fourier space, the velocity fluctuation may be written as, 
5 v (k) — —ffi 2 8 where fi is the cosine of the angle between 
the k and line-of-sight unit vectors (Kaiser 1987), anojf] / = 
dlog 5/dlog (1 + z). To leading order in 5, it then follows 
from equation ^ that the PS of brightness temperature 
fluctuations is given by, 

Pat = T b 2 [(xjaPss - 2x H i(l - x m )Ps x + (1 - x m ) 2 P xx ) 
+2fn 2 (x m Pss - im(l - x m )Psx) + /V (xmPss)] (22) 

where % = 23.8 [(1 + z)/10]' mK, and P ss , P Sx and P xx are 
the PS of density fluctuations, the PS of ionization fluctua- 
tions, and the cross-PS of ionization and density fluctuations 
respectively. Our expression for S x allows the ionization PS 
and ionization-density cross-PS to be related to the density 
PS in a very simple way 



P Sx = 
p 



xm 
1 - xhi 

•r 2 



[b(l-K(k))-l]P ss 
[b{l-K(k))-l] 2 P ss . 



, - ., >>-.. < (23) 

(1 — xm) 

In addition, the ionization PS and ionization-density cross- 
PS are very simply related to each other 



xm 



(1 — xm) 



[b(l-K(k))-l]P Sx 



(24) 



Upon substitution of equation (|23[) into equation (|22l) 
we obtain an expression for the angular dependence of the 
post-reionization 21cm PS 



Pat(6) = T b 2 x| : [6(1 - K(k)) + ffi 2 ] 2 P ss . 



(25) 

Thus, when combined equations (122 p and (|23[) provide a sim- 
ple relation between the density PS and the 21cm PS involv- 
ing ionization terms, allowing the full 21cm PS to be utilized 
in deriving cosmological constraints. This compares favor- 
ably with the reionization epoch, during which either the 
/Li 4 term alone can be utilized (Barkana & Loeb 2005), or a 
functional form for the ionization PS must be assumed (Mao 
et al. 2008; Rhook et al. 2008). On large scales fluctuations 
in the ionization field will be correlated with density fluctu- 
ations. Thus fluctuations in ionizing background would be 
degenerate with the galaxy bias 6. However, on small scales 
there is scale dependence in the function K which breaks the 

8 The quantity / is close to unity at high redshifts, taking values 
of 0.974, 0.988 and 0.997 at z = 2.5, 3.5 and 5.5. 
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degeneracy between ionization and galaxy bias. In addition, 
the angular dependence of the PS, which is due to gravi- 
tational infall and hence is independent of galaxy bias or 
ionization, breaks the degeneracy between galaxy bias and 
neutral fraction. 

Equation (|25|l is valid on scales for which fluctuations 
are in the linear regime, and we can assume that galaxy 
bias is not scale dependent. In the case where b = / = 1 
and K = (i.e. unbiased sources at high redshift, with 
no modulation by the ionizing background) our expression 
for the PS reduces to the form of a uniformly ionized IGM, 
Pat(&) = T b 2 x H [1 + n 2 ] 2 P SS (Barkana & Loeb 2005). We 
note that in the absence of a modulation of power due to the 
ionizing background, equation (|25[) could have been derived 
directly from considering the galaxy PS (once it is realized 
that effect of a line-of-sight velocity is the same for the opti- 
cal depth in the 21cm line and the galaxy number density). 

Equation (I25p is valid for a single halo mass. However 
the DLAs reside within halos having a range of masses. To 
account for this we define an unknown HI mass weighted 
probability density function p(b)db for the galaxy bias. By 
substituting 5 X into equation {l}, and integrating over p(b) 
it is easy to show that the PS for an arbitrary distribution 
of bias is 



Pat = T h 2 x 2 m P M [(b) (1 - K(k)) + ffi 2 ] 



21 2 



(26) 



where (b) is the mean of the bias distribution p(b). For com- 
pleteness we also calculate the spherically averaged PS 



psph 
AT 



= TbXmPss 

{b) 2 (l-K(k)) 2 + ^f(b)(l-K(k)) + ^f 2 



(27) 



which, in the case of (6) = 1, / = 1 and K a = re- 
duces to the form for a uniformly ionized high redshift IGM, 
P A P £ = 1.87 x T^x^jPss (Barkana & Loeb 2005). Note that 
although the relation between neutral hydrogen mass and 
host halo mass is unknown, the PS depends only on the 
mean of p(b). For later use we define M/m to be the halo 
mass corresponding to (b). 

In Figure Q] we show examples of the dimensionless 
21cm PS after reionization [A = k 3 Pat / '(2ir 2 )]. We assume 
Mn,) = lO n M0, use the Sheth-Tormen fitting function for 
the bias (Sheth, Mo & Tormen 2001), and consider three red- 
shifts, z = 2.5, z = 3.5 and z = 5.5, at which the mean-free- 
path for ionizing photons is taken to be A m f p = 600, 300 and 
100 co-moving Mpc, respectively (Bolton & Haehnelt 2007; 
Faucher-Giguere et al. 2008). There is at least a factor of 
2 uncertainty in the value of the mean- free-path. However 
our results are insensitive to the exact value because the 
mean-free-path is larger than the scales accessible to 21cm 
PS observations (owing to the requirements of foreground 
removal) at all redshifts except those closest to reionization. 

Two models are shown to illustrate the effect of the 
ionizing background on 21cm fluctuations. In the first we 
assume a fiducial model for which the ionizing background 
has no effect on the fluctuations in neutral fraction (K — 0). 
This is shown by the dotted curves in Figure [T] We also show 
examples of a model in which Ko = 1, a value which is un- 
physically large but is chosen so as to clearly illustrate the 
qualitative effect of fluctuations in the ionizing background 
on the shape of the 21cm PS. In this model, the fluctuations 
in neutral fraction due to the ionizing background are as 



strong as (and hence cancel) the effect of the galaxy bias on 
large scales. Figure [1] illustrates that the fractional modu- 
lation of the 21cm PS is of order Ko, while as discussed in 
§ 12.2.31 physically plausible values of Ko are of order 10 -2 . 
Since the fractional modification of the PS is of order 2Ko, 
and our results indicate a fractional change in amplitude 
which is smaller than a factor of 2 at k ~ 0.1 over the red- 
shift range considered, we therefore find that fluctuations in 
the ionizing background should modify the 21cm PS after 
reionization by < 1% on the scales and redshifts of interest. 

2.4 Discussion 

The effect of the ionizing background on the 21cm PS can 
be understood qualitatively as follows. Firstly, we have ar- 
gued that the ionizing background follows the galaxy density 
field on large scales. As a result, on large scales regions of 
overdensity (underdensity) will have a higher (lower) than 
average ionizing background. In our linear formulation, the 
effect of the ionizing background is degenerate with galaxy 
bias on large scales, because in the case of Ko > 0, it lowers 
(raises) the contribution to the 21cm intensity from galax- 
ies with a particular halo mass. This modification of the 
21cm intensity of galaxies reduces the amplitude of the PS, 
at all scales k <C fc m f p . We emphasise that because the fluc- 
tuations in the ionizing background are strongly correlated 
with density on large scales, they suppress rather than add 
additional power to the 21cm PS. The suppression becomes 
greater for smaller values of fe/fe m f p . Thus at fixed scale the 
effect of a fluctuating background is larger at higher redshift 
where the mean-free-path is smaller. 

On small scales k 3> fc m f P , only some fraction of the 
ionizing BG is produced locally within the fluctuation. The 
remainder of the ionizing background is generated within a 
region which averages over many fluctuations, and so has a 
mean value that is near that of the background. As a result, 
the suppression of power described above, owing to a cor- 
relation between the ionizing background and overdensity, 
is not as strong on small scales. Our derivation yields the 
variation of the suppression with scale (oc fc _1 ). It should 
be noted that our formulation ignores the possible Poisson 
contribution to fluctuations in the ionizing background due 
to quasars. Poisson fluctuations introduce additional power 
into the 21cm PS beyond the component associated with 
the underlying density field of galaxies. The effect of Poisson 
fluctuations could become important at low redshift (z < 3), 
where quasars contribute significantly to the ionizing back- 
ground. We do not consider the possible Poisson contribu- 
tion in the remainder of this paper which focuses on redshifts 
above the peak of quasar activity (z > 2.5). 

We emphasise that our derivation of the 21cm PS af- 
ter reionization is only valid on scales where the density 
fluctuations are in the linear regime. It should therefore be 
noted that there are potential complications that arise due 
to non-linearities in the mass and velocity fields on small 
scales. Analysis of the galaxy PS derived from N-body sim- 
ulations has shown (Seo & Eisenstein 2005) that the PS 
can be treated as linear on scales greater than 15 co-moving 
Mpc (i.e. fc max <c0.4Mpc -1 ) at z = 3.5, increasing towards 
higher redshifts. However, weak oscillatory features in the 
PS, such as the baryonic acoustic oscillations, are suppressed 
on even larger scales because matter moves across distances 
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Figure 1. Examples of the post-reionization 21cm PS at three redshifts. Two cases are shown, assuming models with K = (dotted 
curves) and K = 1 (solid curves). A DLA host mass of = 1O 11 M0 was assumed at each redshift. The light solid curve shows 

the spherically averaged noise for the MWA5000, assuming lOOOhr integration on each of three fields. The sharp upturn at low k is due 
to the assumption that foreground removal prevents measurement of the PS at scales corresponding to a bandpass larger than 8MHz. 
The thermal noise and cosmic variance (including the Poisson noise) components are plotted as the light dotted and light dashed curves 
respectively. The sensitivity curves are plotted within bins of width Afc = fc/10. 



on the order of ~ 5-10Mpc over a Hubble timc0. As groups 
of galaxies form, the linear-theory prediction for the loca- 
tion of each galaxy becomes uncertain, and as a result noise 
is added to the correlation among galaxies and hence to the 
measurement of the mass PS. The noise associated with the 
movement of galaxies smears out the acoustic peak in the 
correlation function of galaxies in real space (Eisenstein, Seo 
& White 2007; Seo, Siegel, Eisenstein & White 2008). The 
associated reduction of power in the baryonic acoustic os- 
cillations is found to be in excess of 70% on scales smaller 
than fc max ~ 0.4Mpc _ at z = 3.5, corresponding to a length 
scale of ~ 7r/(2fc max ) = 3.9 comoving Mpc (Seo et al. 2008). 

More importantly for experiments that aim to measure 
the PS shape is the non-linear correction to the PS due 
to virialised groups of DLAs (e.g., Tinker et al. 2006; Tin- 
ker 2007). The group velocity dispersion, ovir, introduces the 
so-called "finger-of-god" in redshift space which is unrelated 
to the peculiar velocity according to linear theory. We can 
estimate the wavenumber at which this effect becomes im- 
portant relative to the Hubble flow as, fc v i r ~ it H (z) / (2a v i T ) . 
At z ~ 3.5 we find Aw ~ lOMpc -1 . In the future, an im- 
provement to our analysis could be made by including the 
analytic model for the redshift space galaxy two-point cor- 
relation function described by Tinker (2007). This model 
is constructed within the framework of the Halo Occupa- 
tion Distribution, and quantifies galaxy bias on linear and 
non-linear scales. In addition, the model describes redshift- 
space distortions and clustering on both linear and non- 
linear scales. Finally, other non-linear effects may arise that 
are not present in galaxy redshift surveys, owing to internal 
rotation curve of the neutral hydrogen. However, the corre- 
sponding wave numbers are very large, fc 2> lOMpc -1 . 



9 This characteristic scale of displacement follows from the fact 
that erg, the normalization of the power-spectrum on 8/i — 1 Mpc, 
is of order unity at the present time. 



2.5 sensitivity to the 21cm PS after reionization 

Before proceeding to discuss the cosmological potential of 
the post-reionization 21cm PS, we compute the sensitivity 
with which it could be detected. 

To compute the sensitivity APat(&,m) °f a radio- 
interferometer to the 21cm PS, we follow the procedure out- 
lined by McQuinn et al. (2006) and Bowman, Morales & 
Hewitt (2007) [see also Wyithe, Loeb & Geil (2008)]. The 
important issues are discussed below, but the reader is re- 
ferred to these papers for further details. The uncertainty 
comprises of components due to the thermal noise, and due 
to sample variance within the finite volume of the obser- 
vations. We also include a Poisson component due to the 
finite sampling of each mode (Wyithe 2008), since the post- 
reionization 21cm PS is generated by discrete clumps rather 
than a diffuse IGM. However we find that Poisson noise 
dominates only when M (i)) > lO^M© (see Wyithe 2008). 
We assume that foregrounds can be removed over 8MHz 
bins, within a bandpass of 32MHz (McQuinn et al. 2006) 
[foreground removal therefore imposes a minimum on the 
wave-number accessible of fc ~ 0.07[(1 + z)/4.5] _1 Mpc ], 
and consider a hypothetical follow-up telescope to the MWA 
which would have 10 times the total MWA collecting area 
(we refer to this as the MWA5000). This telescope is as- 
sumed to have an antenna density distributed as p(r) oc r~ 2 
within a diameter of 2km and a flat density core of radius 
80m (see McQuinn et al. 2006). The antennae design is taken 
to be optimized at the redshift of observation (so that the 
physical collecting area of the array equals its effective col- 
lecting area), with each of 5000 phased arrays (tiles) con- 
sisting of 16 cross-dipoles. An important ingredient is the 
angular dependence of the number of modes accessible to 
the array (McQuinn et al. 2006). We assume 3 fields are ob- 
served for 1000 hr each. For the experiment described, the 
signal-to-noise of the PS will be largest over the decade of 
scales around k ~ OTMpc -1 (McQuinn et al. 2006; Wyithe, 
Loeb & Geil 2008). The sensitivity curves (within bins of 
Afc = fc/10) are plotted as the solid grey lines in Figure [1] 
We find that the both the 21cm PS, as well as the effect of 
a non-zero value of K would be easily measured using the 
MWA5000. 
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We note that the unprecedented power of cosmic- 
variance limited 21cm surveys for constraining cosmologi- 
cal parameters is made possible by the wide fields of view 
for the low-frequency telescopes under construction, com- 
bined with the very large volumes available at high red- 
shift (Loeb & Wyithe 2008; Mao et al. 2008). For ex- 
ample, in units of the Sloan Digital Sky Survey volume 
(Vaioan ~ 5.8 x 10 8 Mpc 3 ) the telescope described would ob- 
serve V/V s i oan ~ 4.8, V/V s i oan ~ 7.6 and V/V s i oan ~ 9.8 per 
observing field, at z = 2.5, z = 3.5 and z = 5.5, respec- 
tively. These large volumes are obtained because of the very 
wide field of view available to an array with a design like 
the MWA. The array described would have a primary beam 
of A$ ~ 2000 square degrees. 



3 APPLICATION OF THE 

ALCOCK-PACZYNSKI TEST 

Like traditional galaxy redshift surveys, the observed 21cm 
PS is sensitive to only one underlying PS (Pss), in addition 
to the 4 parameters, xm and (6), fc m f p and K . These param- 
eters are related to properties of DLAs and LLSs, and can be 
measured through independent means via quasar absorption 
line studies (particularly true for xva, (b) and fc m f p ). This sit- 
uation should be contrasted with the reionization era where 
the observed 21cm PS is also sensitive to long range, non- 
gravitational fluctuations through Pg a and P xx . Provided 
that the non-linear evolution of the PS can be properly ac- 
counted for, this dominant dependence on the density PS 
provides a considerable advantage for the purposes of cos- 
mological constraints. 

To illustrate the potential for cosmological constraints 
from the anisotropy of the 21cm PS after reionization, wc 
calculate the Alcock-Paczynski (AP) effect. Our approach is 
to specify the general result of Barkana (2006) for the distor- 
tion of the true PS [Pat(^>aO] resulting from an incorrect 
choice of cosmology. This is parameterised in terms of the 
dilation parameters a and a±, which describe the distor- 
tions between the transverse and line-of-sight scales, and in 
the overall scale, respectively. These are defined such that 
(1 + q) is the ratio between the assumed and true values 
of (DaH), while (1 + a±) is the ratio between the assumed 
and true values of the angular diameter distance, Da- In 
the AP test, the correct cosmology is inferred by finding 
cosmological parameters for which a — a± = 0. 

To calculate the AP effect we apply equation (8) in 
Barkana (2006) 



P&T{k,n) 



(1 + a — 3a±)P A T + (otfi 
2\ dP\ T 



dP. 



AT 



dink 



+a(l — fi 



9 In /i 



(28) 



to the PS in equation (|26|l . Here the PS Pat and Pat are 
evaluated at the observed k. This procedure results in a 



modified PS that is related to the true density PS (Pss) via 
Pat(*,p) = T h 2 x 2 m Ph((b)(l-K(k)) + ffi 2 ) 2 

dink 

T 2 x 2 lI P t ss ((b)(l-K(k)) + fn 2 ) 



-K(k) 



■(b) (afi 



1 + (k/k m tp) 

- 4a (1 - n 2 ) ffi 2 ] ■ (29) 

Barkana (2006) expanded the general form of this equation 
to find the linear combinations of the true PS Pss, Ps x an d 
P xx which form the coefficients of terms containing fi°, fi 2 , 
fj, 4 and n 6 . His result shows that in general, the parameter a 
(which corresponds to anisotropy) mixes Pg S , Pg x and P xx . 
As a consequence, the fj. 6 term, which arises through the 
AP effect, must be isolated in order to probe the anisotropy 
parameter a. On the other hand, equation (|29f) is sensitive 
only to Pss- As a result, the coefficients of all powers of jj, in 
the post-reionization PS can be utilised in the AP effect to 
find the cosmological parameters that yield a = aj_ = 0. 



3.1 AP constraints on the PS 

We next use equation (|29jl to calculate the permissible region 
of parameter space p = (a, a±, xm, Mm , K ) around a true 
solution with PS P% s and p a = (0, 0, 0.02, 10 U M Q , 0). The 
fiducial value of Mm is chosen to lie in the middle of the 
mass range for DLAs measured from the cross-correlation 
analysis of Cook et al. (2006). We have assumed that the 
mean-free-path of ionizing photons is known a'priori, and 
do not fit it as a free parameter. As part of this procedure, 
we assume an estimate for the sensitivity of the future low- 
frequency interferometer described in § 12.51 to the 21cm PS 
[APat(&, aO]i and construct likelihoods 

^--sE ^^ffi**' )'- (30) 

where the sum is over bins of k and fi. To estimate the 
constraints achievable via the AP effect, we modify the lin- 
ear power spectrum to include the nonlinear erasure of the 
acoustic peaks based on the work of Eisenstein et al. (2006) . 
Specifically, we use 



Pss(k) = (Pss(k) 



pir h (k) 



exp(-fc 2 E? ll /2)+P|r h (fc), 



(31) 

where P|™ th is the "no wiggle" form from Eisenstein & Hu 
(1999), and E ni = 3.9[(1 + z)/4.5] _1 Mpc in the high red- 
shift limit. In addition, to account for the possibility of non- 
linearity in the smooth PS at small scales we restrict our 
fitting to wave numbers fc max < O^Mpc -1 (see § I2.4[I , 

The results are shown in Figure [2] for three redshifts, 
z — 2.5, 3.5 and 5.5. Likelihood contours are shown at 
each redshift for the parameter sets (a, a± ) , (xm, Mm ) and 
(K ,M(h))- For each set the likelihood is marginalised over 
the remaining 3 parameters assuming flat prior probabili- 
ties. The exception is Shi ~ 0.02 ± 0.002, which we have 
assumed to be known a'priori (with Gaussian errors) from 
future observations of quasar absorption spectra. The as- 
sumed uncertainty in xhi corresponds to a factor of ~ 2 
improvement over existing measurements (Wolfe, Gawiser 
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Figure 2. Constraints on the post-reionization PS distortions achievable via the AP test at 3 redshifts, assuming a fiducial model with 
K = 0. The three columns show parameter sets (a, a±), (xm, M^)) an d (Ko, ^(b))- I n each case likelihood contours are shown at 60% 
and 7% of the maximum likelihood. The projections of these contours onto individual parameter axes correspond to the 68.3% and 90% 
confidence ranges respectively. 



& Prochaska 2005). Figure [2] shows that the angular de- 
pendence of the PS breaks degeneracy between the different 
parameters, allowing xhi, , a, and aj_ to each be mea- 
sured from the fitting. In addition, the departure from the 
underlying shape of the PS allows the value of K D to also be 
constrained. 

We note that by including the uncertainties in the AP 
effect (through a and a± ) , our derived uncertainties for non- 
cosmological parameters (M/fA, xhi and K D ) include the un- 
certainty in the underlying matter power spectrum. An ex- 
ception is that the uncertainty in the PS amplitude (which 
is proportional to the normalization of the primordial PS, 
as) is degenerate with zhi. However, since the error in xhi 
is comparable to or larger than the fractional error in erg, 
this does not add additional uncertainty relative to the con- 
straints shown in Figure [2] An important additional uncer- 
tainty may also be introduced through the uncertainty in 
the primordial PS index or through a running spectral in- 
dex, which we have not considered in this study. 

At 2 < 3.5 the neutral fraction constraints are improved 
over the assumed prior information from the IGM. The mass 
of the DLA hosts is constrained to high precision (tens 



of percent, Wyithe 2008), with more accurate estimates at 
lower redshifts. In addition, the value of K a is constrained 
to be smaller than a few hundredths (with minimal redshift 
dependence), although this in not at a level comparable to 
the physically expected value of ~ 10~ 2 . Thus, on scales 
larger than the mean-free-path, the perturbation due to the 
ionizing background ((b) K S) is constrained at a level that 
is smaller than (S 2 ) 1 ^ 2 . On scales smaller than the mean 
free path, this perturbation is suppressed by (fe/fe m f p ) . 
The cosmological constraints follow from the precision with 
which the distortion of the 21cm PS (as described by a and 
a±) can be measured. We find that distortions owing to the 
assumption of an incorrect cosmology could each be con- 
strained at a percent level. 



3.2 Comparison with high redshift constraints 

Before concluding we show for comparison (Figure [3} the 
constraints on the parameter set (a, ax) for optimistic 
high redshift cases (z = 12.5 and 15), where we assume 
thi = 1 and fluctuations dominated by the density field 
(Pxx = Psx = 0). These cases have fitted parameters 
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Figure 3. Central and Right panels: Constraints on the pre-reionization PS distortions achievable via the AP test at two redshifts. The 
results correspond to an optimistic high-rcdshift case, with Sjji = 1; an d Pxx = Pgx = 0- The panels show the parameter set (a, ctj_). In 
each case, likelihood contours are shown at 60% and 7% of the maximum likelihood with each set marginalised over Sjji (which we have 
assumed to be unconstrained a'priori). The projections of these contours onto individual parameter axes correspond to the 68.3% and 
90% confidence ranges respectively. For comparison, (left panel) the K = case at z = 3.5 is repeated from Figure[2] 



p— (a, a±,xm), and correspond to the IGM prior to the 
formation of galaxies and ionized bubbles. The smaller xhi 
at z < 6 relative to the pre-reionization epoch is offset by 
the galaxy bias factor enhancement as well as the growth of 
structure, and by lower foreground contamination. With re- 
spect to constraining cosmological parameters, we find that 
the post-reionization 21cm PS would be competitive with 
the most optimistic expectations for the 21cm PS at high 
redshift. 



4 CONCLUSIONS 

In this paper we have derived the 21cm power spectrum 
(PS) following the completion of reionization. Our approach 
is to derive the PS within the formalism that has been devel- 
oped to calculate the 21cm PS during reionization (Barkana 
& Loeb 2005; McQuinn et al. 2006; Mao et al. 2008). By 
framing the derivation this way we are able to directly com- 
pare the relative merits of cosmological constraints from 
21cm surveys at redshifts prior to and post reionization. We 
have derived expressions for the components of the post- 
reionization 21cm PS that are due to the density field, the 
ionization field, and their cross-correlation. As is the case 
prior to reionization, we find that these components con- 
tribute to the observed PS in proportions that depend on 
the angle relative to the line-of-sight along which the power 
is measured. However, in difference from the situation prior 
to reionization, we have shown that all components of the 
21cm PS are directly proportional to the PS of the un- 
derlying matter fluctuations, with a small but predictable 
modulation on scales below the mean-free-path of ionizing 
photons. We have derived the form of this modulation, and 
have shown that its effect on the observed PS will be at 
less than the 1% level for physically plausible astrophysi- 
cal models of DLA systems. The simplicity of the 21cm PS 
after reionization stands in contrast to the astrophysical un- 
certainty during the reionization epoch, where HII regions 
dominate the 21cm signal. This simplicity will enhance the 
utility of the 21cm PS after reionization as a cosmological 
probe (Loeb & Wyithe 2008) by removing the need to sep- 
arate the PS into physical and astrophysical components 
(Barkana & Loeb 2005). 



To illustrate the utility of the 21cm PS after reioniza- 
tion as a cosmological probe, we have examined the Alcock- 
Paczynski (1979) test. Our calculations show that the next 
generation of low-frequency arrays could measure the angu- 
lar distortion of the PS to around ~ 1% at z ~ 3.5. It has 
previously been shown that the scale of Baryonic Acous- 
tic Oscillations, which constitutes a standard ruler (Blake 
& Glazebrook 2003; Seo & Eisenstein 2005; Eisenstein et 
al. 2005; Padmanabhan et al. 2007), can be also be used to 
independently probe H and Da in 21cm PS, and hence to 
measure the equation of state of the dark energy (Wyithe, 
Loeb & Geil 2008; Chang et al. 2007). With the caveat that 
non-linear evolution of the 21cm PS must be quantitatively 
understood, our simple analysis indicates that the precision 
achievable via the Alcock-Paczynski (1979) test using the 
21cm PS after reionization could be better than those avail- 
able via a 21cm fluctuation measurement of the acoustic 
scale of baryonic oscillations for a given observing strategy 
(Shoji, Jeong & Komatsu 2009). More generally, our anal- 
ysis shows that the 21cm PS after reionization shares the 
same favorable features as a galaxy redshift survey. The ad- 
vantage of using 21cm fluctuations lies in the fact that in- 
dividual sources need not be resolved. This would allow PS 
measurements using 21cm fluctuations to be extended to 
higher redshifts. 

Acknowledgments. The research was supported by the 
Australian Research Council (JSBW), by NASA grants 
NNX08AL43G and LA, by FQXi, and by Harvard University 
funds (AL). We thank an anonymous referee for correcting 
an earlier error. 

REFERENCES 

Alcock C. & Paczynski, B., 1979, Nature, 281, 358 
Barkana, R., 2006, Mon. Not. R. Astr. Soc, 372, 259 
Barkana, R. & Loeb, A., 2007, Reports of Progress in 
Physics 70, 627. 

Barkana, R. & Loeb, A., 2005, Astrophys. J. Lett. 624, L65 
Blake C. & Glazebrook, K., 2003, Astrophys. J. 594, 665 
Bolton, J. S., & Haehnelt, M. G. 2007, MNRAS, 382, 325 
Bowman, J.D., Morales, M.F. & Hewitt, J.N., 2007, Astro- 
phys. J., 661, 1 



© 0000 RAS, MNRAS 000, 000-000 



10 Wyithe & Loeb 



Chang, T.-C, Pen, U.-L., Peterson, J.B. & McDonald, P., 

2007, ArXive-prints, 709, larXiv:0709.3672l 

Cooke, J., Wolfe, A. M., Gawiser, E., & Prochaska, J. X. 
2006, ApJ, 652, 994 

Corbelli, E., Salpeter, E. E., & Bandiera, R. 2001, ApJ, 
550, 26 

Curran, S.J., Tzanavaris, P., Pihlstrom, Y.M. & Webb, 
J.K., 2007, Mon. Not. R. Astron. Soc., 382, 1331 
Eisenstein, D. J. et al., 2005, Astrophys. J. 633, 560 
Eisenstein, D. J., & Hu, W. 1999, ApJ, 511, 5 
Eisenstein, D. J., Seo, H.-J., & White, M. 2007, ApJ, 664, 
660 

Faucher-Giguere, C, Lidz, A., Hernquist, L., & Zaldar- 
riaga, M. 2008, ArXiv e-prints, 807, larXiv:0807.4T77l 
Furlanetto, S. R., Oh, S. P. & Briggs, F. H., 2006, Phys. 
Rep. 433, 181. 

Iliev, I.T., Mellema, G., Pen, U.-L. & Shapiro, P.R., 2007, 

ArXiv e-prints, 712, larXiv:0712. 13561 

Kaiser, N., 1987, Mon. Not. R. Astron. Soc, 227, 1 

Komatsu, E., et al. 2008, ArXiv e-prints, 803, 

larXiv:0803.0547l 

Lewis, A. & Challinor, A., 2007, Phys. Rev. D 76, 083005 
Loeb, A. & Zaldarriaga, M., 2004, Phys. Rev. Lett. 92, 
211301 

Loeb, A. & Wyithe, J.S.B., 2008 Physical Review Letters, 
100, 161301 

Mao, Y., Tegmark, M., McQuinn, M., Zaldarriaga, M. & 
Zahn, O., 2008 ArXiv e-prints, 802. larXiv:0802.1710l 
McQuinn, M., Zahn, O., Zaldarriaga, M., Hernquist, L. & 
Furlanetto, S.R., 2006, Astrophys. J. 653, 815 
Osterbrock, D. E., & Ferland, G. J. 2006, Astrophysics 
of gaseous nebulae and active galactic nuclei, 2nd. ed. by 
D.E. Osterbrock and G.J. Ferland. Sausalito, CA: Univer- 
sity Science Books, 2006, 

Padmanabhan, N. et al., 2007, Mon. Not. R. Astr. Soc, 
378, 852 

Pen, U.-L., Chang, T.-C, Peterson, J. B., Roy, J., Gupta, 
Y., & Bandura, K. 2008, The Evolution of Galaxies 
Through the Neutral Hydrogen Window, 1035, 75 
Pen, U.-L., Staveley-Smith, L., Peterson, J. & Chang, T.- 
C, 2008, ArXiv e-prints, 802. larXiv:0802.3239l 
Pritchard, J.R. & Loeb, A., 2008, ArXiv e-prints, 802, 
larXiv:0802.2T02l 

Rhook, K.J., Geil, P.M. & Wyithe, J.S.B., 2008, ArXiv e- 
prints, 803, larXiv:0803.2802l 

Santos, M.G., Amblard, A., Pritchard, J., Trac, H., Cen, R. 
& Cooray, A., 2007, ArXiv e-prints, 708. larXiv:0708. 24241 
Schaye, J. 2004, ApJ, 609, 667 

Shoji, M., Jeong, D., & Komatsu, E. 2009, ApJ, 693, 1404 
Seo, H.-J. & Eisenstein, D.J., 2005, Astrophys. J. 633, 575 
Seo, H.-J., Siegel, E. R., Eisenstein, D. J., & White, M. 

2008, ArXiv e-prints, 805, larXiv:0805.0TT7l 

Shoji, M., Jeong, D., & Komatsu, E. 2008, ArXiv e-prints, 
805. larXiv:0805.4238l 

Sheth, R.K., Mo, H.J. & Tormen, G., 2001, Mon. Not. R. 

Astron. Soc, 323, 1 

Susa, H. 2008, ApJ, 684, 226 

Tinker, J. L. 2007, MNRAS, 374, 477 

Tinker, J. L., Weinberg, D. H., & Zheng, Z. 2006, MNRAS, 
368, 85 

Wolfe, A.M., Gawiser, E., & Prochaska, J.X., 2005, Ann. 
Rev. Astr. Astrophys. 43, 861 



Wyithe, J.S.B. & Loeb, A., 2008 Mon. Not. R. Astr. Soc, 
383, 606 

Wyithe, J.S.B., Loeb, A. & Geil, P.M., 2008, Mon. Not. R. 

Astr. Soc, 383, 1195 

Wyithe, S. 2008, ArXiv e-prints, 804, larXiv:0804. 16241 
Zheng, Z., & Miralda-Escude, J. 2002, ApJL, 568, L71 



© 0000 RAS, MNRAS 000, 000-000 



